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Abstract 

Context. Shock fronts in young supernova remnants are the best candidates for being sites of cosmic ray acceleration up to a few PeV, 
though conclusive experimental evidence is still lacking. 

Aims. Hadron acceleration is expected to increase the shock compression ratio, providing higher postshock densities, but X-ray 
emission from shocked ambient medium has not firmly been detected yet in remnants where particle acceleration is at work. We 
exploited the deep observations of the XMM-Newton Large Program on SN 1006 to verify this prediction. 

Methods. We performed spatially resolved spectral analysis of a set of regions covering the southeastern rim of SN 1006. We studied 
the spatial distribution of the thermodynamic properties of the ambient medium and carefully verified the robustness of the result with 
respect to the analysis method. 

Results. We detected the contribution of the shocked ambient medium. We also found that the postshock density of the interstellar 
medium significantly increases in regions where particle acceleration is efficient. Under the assumption of uniform preshock density, 
we found that the shock compression ratio reaches a value of ~ 6 in regions near the nonthermal limbs. 

Conclusions. Our results support the predictions of shock modification theory and indicate that effects of acceleration of cosmic ray 
hadrons on the postshock plasma can be observed in supernova remnants. 

Key words. X-rays: ISM - ISM: supernova remnants - ISM: individual object: SN 1006 



1. Introduction 

Supernova remnant (SNRs) are strong candidates for being the 
main source of energetic cosmic rays up to at least 3 x 10 15 eV 
dBlandford & Eichlerfl987l:lBerezhko & Volkl2007l) . X-ray syn- 
chrotron emission from high-energy electrons accelerated at the 
shock front up to Te V energies was first observed in SN 1006 
( Kova ma et a fl [l99l and then in other young SNRs (Reynolds 
2008), but a firm detection of high-energy hadrons is still lack- 
ing. 

has 



been detected in SN 1006 
handful o f SNRs show- 
e mission ( Aharonian et al J 
lAharonian etail l2009t 
H.E.S.S. Collaboration et alj|201 ll) . The origin of the gamma- 



Recently , TeV emission 
dAcero et all [2010) and in a 
ing b right nontherma l X-ray 
200 7albT: lAlbert et all 1200' 



ray emission can be leptonic (i. e., inverse Compton from the 
accelerated electrons) and/or hadronic (i. e., proton-proton 
interactions with n production and subsequent decay). The 
hadronic scenario would directly prove that SNRs can accel- 
erate cosmic rays up to PeV energies. Unfortunately, it is not 
easy to unambiguously ascertain the origin of the gamma-ray 
emission, as, for example, in RX J17 13.7-3946, where both 
hadro nic dBerezhkp & Volkl l2010h and lepto nic dAbdo et alJ 
l20Tll lEllison et alJl2010h iPatnaude et alJ l2010) scenarios have 
been invoked to explain the observed multiwavelength observa- 
tions. However, different indications concur in supporting the 
presen ce of high-energy hadrons in some young SN Rs (e. g., 
Tvcho. lEriksen et al.ll201 lHMorlino & Capriolill2012l and RCW 



86. lHelder et~ai1l2009l) . As for SN 1006, a pure leptonic model 
is consistent with the observations, but a mixed scenario that 
includes leptonic and hadr onic components a lso provides a good 
fit to the gamma-ray data dAcero et a l. 2010). 

An alternative way to reveal hadron acceleration in SNRs 
is to probe its effects on the shock dynamics. The non- 
linear back-reaction of high-energy particles on background 
plasma is predicted to strongly modify the shock properties 
by increasing the shock compression ratio and decreasing the 
postshock tempe rature with respect to the expected Rankine- 
Hugoniot values (|Berezh ko & Ellison 1999: iDecourchelle et al.l 
2000; iBlasil 120021: IVink et alJl201dh . This effect is known as 
"shock modification". The observational test of these predic- 
tions requires accurate diagnostics of the thermal X-ray emis- 
sion from the shocked interstellar medium (ISM). However, 
the ISM contribution in the X-ray spectra of remnants where 
particle acceleration is efficient is typically masked out by the 
bright synchrotron emission and by th e thermal emission from 
sh ocked ejecta, as in R XJ17 13.7-3 946 dAcero et al.ll2009t) Vela 
Jr dPannuti et alJl2010h. Gl. 9+0.3 dBorkowski et alj feoiO). and 
Tvcho dCassam -Chenai et al. 2007), and has not been firmly de- 
tected so far. In SN 1006 the quest for X-ray emission from 
shock ed ISM is also ongoing dAcero et al.l 120071: iMiceli et al.l 
120091) . 

SN 1006 exhibits a morphology characterized by two op- 
posed radio, X-ray, (and 7— ray) bright limbs dominated by 
nonthermal emission (bilateral SNR) and separated by an in- 
ner region of low surface brightness and soft, thermal X-ray 
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emission. Ther mal X-ray emission has been asso ciated with 
shocked ejecta (lAcero et al.l l2007: Mice li et alJ20 09). consistent 
with the detection of F e-rich plasma in the interior of the shell 
dYamaguchi et alJ2008l) . A model with a nonthermal component 
plus a thermal component associated with the ejecta can fit the 
archive XMM-Newton spectra. An additional thermal ISM com- 
ponent is not needed from a statistical point of view (lAcero et alJ 
120071: iMiceli et al.ll2009h . though some hints about ISM temper- 
ature and d ensity (£7\sm ~ 1-5 - 2 keV, n < 0.2 cm 4 ) can 
be inferred (lAcero et alJl2007h . Suzaku spectra of SN 1006 have 
been modeled with three thermal components only by assuming 
that the soft component originated in the ISM; nevertheless, it 
was not possible to exclude that the O lines, which dominate the 
soft emission, originate in the ejecta ( Yamaguchi et al. 20Q8|). A 
firm detection of X-ray emission from the ISM is therefore still 
lacking. 

Indirect evidence for shock modification in SN 1006 has 
been obtained by measuring the distance, D$fcd, between the 
shock front and the contact discontinuity, which is expected to 
be smaller in nonthermal limbs where particle acceleration is 
more efficient. D$fcd is instead almost the same all over the 
shell (even in regions dominated by thermal emission), though 
it is much smaller than that expected from unm odified shocks 
dCassam-Chenai et all 120081: IMiceli et alj|2009h . Recently, 3-D 
magneto-hydrodynamic simulations have shown that this small 
distance can be naturally expl ained by ejecta clum ping, without 
invoking shock modification dOrlando et al.l |20T2I) . Therefore, 
the small value of D$ fcd is not a reliable indicator of hadronic 
acceleration. 

Here new, deep observations of the XMM-Newton SN 1006 
Large Program allow us to present the detection of X-ray emis- 
sion from shocked ISM in SN 1006 and to show indications for 
the effects of hadron acceleration on the postshock density of the 
ambient medium. 

The paper is organized as follows: in Sect.|2l we describe the 
data analysis procedure and the background subtraction; in Sect. 
[3] we show the results of the spatially resolved spectral analysis; 
and, finally, we discuss our conclusions in Sect. [4] 

2. Data analysis and background subtraction 

In this work, we analyze the data obtained within the XMM- 
Newton Large Program of observations of SN 1006 (PI A. 
Decourchelle, 700 ks of total exposure time). These observa- 
tions were all performed with the medium filter by using the full 
frame mode for the MOS cameras and the extended full frame 
mode for the pn camera, and were processed with the Science 
Analysis System. Light curves, images, and spectra were cre- 
ated by selecting events with PATTERN<12 for the MOS cam- 
eras, PATTERN <4 for the pn camera, and FLAG=0 for both. 
To reduce the contamination by soft proton flares, the origi- 
nal event files were screened by using the sigma-clipping algo- 
rithm (ESPFILT tasks). The ima ges were pro d uced by adopt- 
ing the procedure described in Mic eli et alJ d2009l) and are 
all background-subtracted, vignetting-corrected, and adaptively 
smoothed to a signal-to-noise ratio of ten. 

Spectral analysis was performed in the energy band 0.3 - 7 
keV using XSPEC. The source regions a - g (see Sect. [3] and 
Fig. Q} analyzed in this paper are covered by the observations 
ID 0555630101 and 0555631001, while region h is covered also 
by observation ID 0653860101 . The screened MOS 1 /MOS2/pn 
exposure times are 43/42/26 ks, 59/60/43 ks, and 101/104/71 
ks for observations 0555630101, 0555631001, and 0653860101, 
respectively. MOS and pn spectra of the different observations 



were fitted simultaneously. Pn spectra were analyzed only in the 
0.6 - 3 keV band to reduce possible contamination from the in- 
strumental background (that can be high above 3 keV) and cross- 
calibration issues (below 0.6 keV). The results of the fittings of 
the pn spectra in the 0.6-7 keV band are shown in the Appendix. 

The spe c tral m odel adopted here is the same as that in 
IMiceli et al] d2009l) (VPSHOCK component for the ejecta plus 
nonthermal emission) with the addition of a VPSHOCK com- 
ponent with solar abundances to describe the ISM. VPSHOCK 
model dBorkowski et aljfeOOll) describes an isothermal and op- 
tically thin plasma in nonequilibrium of ionization with a linear 
distribution of ionization timescale versus emission measure (its 
parameters are electron temperature, plasma emission measure, 
and ionization parameter). The nonthermal e mission is described 
by the SRCUT dRevnolds & Keohand 19991) model (synchrotron 
emission from an electron power law with exponential cut-off). 
In the SRCUT c omponent, we fixed the radio-to-X- ray photon 
index to a = 0.5 dAllen et all2008LlMiceh et al. 2009) , while the 
norm alization was derived from the radio image ( Micel i et al.l 
2009). We verified that the results do not change significantly 
by letting a free to vary in the fittings, though error bars grow 
higher. In the ejecta component, the S abundance is the same as 
the Si abundance. We included three narrow Gaussians in our 
model to take account of the Ks, K e , and K% O VII lines that are 
not included in the spectral code. We also added a systematic 
2% error term to reflect the estimated uncertainties in the cali- 
bration of the instrumental effective area (see the EPIC status of 
calibration report XMM-SOC-CAL-TN-0018). The interstellar 
absorption is described by the TBABS model, and the absorb- 
ing column is set t o Nh — 7 x 10 20 crrT 2 , in agreement with 
iDubner et al. ( 2002). We verified that the best-fit values (and par- 
ticularly the ISM emission measure) do not change significantly 
by letting the Nh vary freely in the fitting procedure (though the 
uncertainties in the best-fit parameters become higher). 

To check the validity and the robustness of our result, we 
adopted three different procedures for background subtraction: 

1 . The best-fit values presented in Sect. [3]and discussed in Sect. 
3] were obtained by subtracting a background spectrum ex- 
tracted from a nearby region immediately outside of the SNR 
shell, BG tot , from the spectrum extracted from our source re- 
gion, SP (source + background). By extracting BG, , from 
different regions outside of the shell, we verified that the 
best-fit values do not depend on the choice of the background 
region since they do not vary significantly. 

2. We also adopted a "double background-subtraction" tech- 
nique: we considered the instrumental background, BG, (ex- 
tracted from the same region in detector coordinates as SP) 
in filter- wheel closed data and computed SP - (BG tot - 
B BGi) - B BGj, where B is the ratio of the flux from the 
unexposed corners of archived observations to that of our 
SN 1006 observations. 

3. We produced the response files for the BG tot spectrum. We 
then modeled the "background" spectra in XSPEC and fixed 
these background parameters when fitting the SP spectra. 

We verified that all the procedures provide consistent results. 
In particular, the values of the ISM temperature obtained with 
the different methods are all consistent within two sigmas, while 
those of the ISM emission measure are all consistent at less than 
one sigmcQ. 



1 The different subtraction methods slightly affect the high-energy 
tail of the spectrum, and small variations of the signal at high energy 
are modeled in the fittings with a higher/lower temperature of the ISM 



2 



M. Miceli et al.: XMM-Newton evidence of shocked ISM in SN 1006: indications of hadronic acceleration 



We used the best-fit values of temperature and emission mea- 
sure to derive the ISM and ejecta densities shown in Sect.[3]and 
th e filling factor (/ism, f eiecta) of each component, obtained as 
in Bocchino et al. ( 1999). These quantities were derived by as- 
suming a pressure equilibrium between ISM and ejecta and by 
calculating the volume of the X-ray emitting plasma (according 
to the formula presented in the Appendix). 



3. Results 

The bilateral morphology of SN 1006 suggests regions at the 
rim with highly efficient particle acceleration (the nonthermal 
limbs), separated by regions with less efficient particle accelera- 
tion. We therefore expect to observe stronger shock modification 
near the nonthermal limbs than in the thermal regions. 

FigureQ]shows an image of SN 1006 in the 2-4.5 keV band 
(upper panel), where the X-ray emission is mainly nonthermal, 
and in the 0.5 - 0.8 keV band (lower panel), where thermal emis- 
sion dominates. The figures also shows the symmetry axis of the 
remnant and the regions selected for the spatially resolved spec- 
tral analysis. The regions cover the southeastern thermal limb 
of SN 1006, where the preshock density of the ISM is expected 
to be fairly uniform (as indicated by the circular sha pe of the 
shock front and in agreement with lDubner et al.ll2002l) . We first 
focus on the middle of the southeastern thermal rim (region e in 
Fig.Q}, where the contribution of the synchrotron emission is the 
smallest and where we expect minimum shock modification. 

3.1. Detection of the ISM component 

Using the deep observations of the Large Program, we found that 
a model with three components (ejecta, ISM, and synchrotron 
emission) fits the spectrum of region e significantly better than 
a model with only two components. A thermal component for 
the ejecta (optically thin plasma in nonequilibrium of ionization 
with free O, Ne, Mg, and Si abundances, plus the nonthermal 
SRCUT component, yields a x 1 = 1 199.8 (with 658 d. o .f.). By 
adding another VPSHOCK component with solar abundances, 
we get a lower value (x 1 = 1112.0 with 655 d. o .f.). The ad- 
ditional component is significant at > 1 1 sigmas. We carefully 
verified the robustness of this result with respect to the analysis 
method and the background subtraction procedure, as explained 
in Sect. [2] The upper panel of Fig. [2] shows the spectrum of re- 
gion e with its best-fit model. 

The origin of the additional thermal component can be at- 
tributed to shocked ISM. This is supported by the values of the 
parameters derived from the spectral analysis: i) the postshock 
density «ism = 0.14 ± 0.01 cm -3 indicates a preshock density 
«o < 0.05 cm" 3 , consistent with the high galactic latitude of the 
remnant and in a greement wit h the upper limit derived from pre- 
vious X-ray data dAcero et alj2007l) : ii) the ionization parameter 
(i. e., the integral of density over time since the passage of the 
shock front) tism ~ 7 x 10 8 s cm 4 provides a very reasonable 
estimate of the time elapsed after the shock impact (~ 200 yr), 
considering the 1 kyr age of the remnant; and iii) the electron 
temperature (~ 1.4 keV) suits the expectations for the shocked 
ISM, as shown below. In collisionless shocks at high Mach num- 
ber (shock speed of the order of 10 3 km/s), the e lectron-to-proton 
temperature ratio is exp ected to be T e /T p <K 1 (Ghavamia n et al.1 
120071: IVink et al.ll2003h . In the northwestern limb of SN 1006, 
the shock velocity is v s f, ~ 3000 km/ s and it has been calculated 




component (and a higher/lower cut-off frequency for the SRCUT com- 
ponent). 



Figure 1. Upper panel: mosaicked count-rate images (MOS- 
equivalent counts per second per bin) of SN 1006 in the 2 - 4.5 
keV band. The bin size is 4" and the image is adaptively 
smoothed to a signal-to-noise ratio ten. The regions selected for 
the spectral analysis of the rim are superimposed (four in the up- 
per panel and four in the lower panel for clarity). North is up and 
East is to the left. The two dashed lines indicate the symmetry 
axis of the remnant, marking the center of the synchrotron limbs 
(northeast and southwest) and of thermal limbs (northwest and 
southeast). Lower panel: same as upper panel in the 0.5 - 0.8 
keV energy band. 



that T e /Tp < 0.07 dGhavamian et al.ll2002l) . Assuming the same 
ratio (and no or little shock modification) in th e southeastern 
regio n, where the shock speed is ~ 5000 km/s dKatsuda et al.l 
2009), we get kT e < 3.5 keV, in agreement with the tempera- 
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channel energy (keV) 



_ SRCUT 



Region "a" 




0.5 1 2 5 

channel energy (keV) 

Figure 2. Upper panel: Spectra of region e (see Fig. [TJ with the 
corresponding best-fit model (see text). The contribution of each 
thermal (ISM, solid line, and ejecta, dashed line) and nonthermal 
(SRCUT, dotted line) components is shown. Only MOS spectra 
are shown for clarity. Lower panel: same as upper panel for re- 
gion a. 



ture above. Moreover, a model developed for RXJ1713. 7-3946 
with parameters similar to those of SN 1006 (but different mass 
of ejecta) and including the effects of particle acceleration and 
electron heating by Coulomb collisions with shocked protons 
has shown that kT e is expected to ra pidly reach 1 keV before 
leveling out between one and two ke V ( Ellis on et al .1201 Oh . in re- 
markably good agreement with our results. Similar values of the 
electron-to-proton temper ature ratio have also been o btained by 
modeling the Tycho SNR dMorlino & Capriolill20"l2h . We there- 
fore conclude that the additional thermal component originates 
in the shocked ISM. 

The upper panel of Fig. [3] shows the confidence contours of 
the ISM density and the synchrotron cut-off frequency obtained 
from the spectrum of region e. In contrast to previous observa- 
tions dAcero et al.ll2007l iMiceli et al. 2009), it is now possible to 
rule out the scenario where the ISM component is not present 
and the cut-off frequency is high (ybreak > 3 x 10 15 Hz). The 
lower panel of Fig. [3] shows the confidence contours of the ISM 
density and the ISM ionization parameter in region e. Since both 
these values depend on the plasma emission measure, it is pos- 
sible to estimate the time elapsed after the shock impact. The 



figure includes isochrones in the («ism, tism) space at 100 yr and 
250 yr. 

3.2. Spatial distribution of the plasma properties 

If efficient hadron acceleration is at work, we expect to observe a 
higher shock compression ratio and a lower shock temperature in 
regions closer to the nonthermal limbs. From the modeling of X- 
ray spectra, we can derive the electron temperature. Howe ver, it 
is stil l unclear how to relate T e with the shock temperature dVinkl 
2012). We therefore consider as the most reliable indicators of 
hadron acceleration the variations in the ISM postshock density 
that can be directly derived through our spectral analysis (see 
Sect. |2j. To quantify the possible shock modification along the 
rim, we analyzed the spectra extracted from the partially over- 
lapping regions [a, b, c, d], and [/, g, h] (shown in Fig.[TJ that 
are closer than e to the nonthermal limbs. 

We first verified that we do not get a significant reduction of 
the^- 2 by leaving the chemical abundances in regions a, b, c, d 
as free parameters in the fitting. Therefore, we fixed the O, Ne, 
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Figure 3. Upper panel: 68%, 90%, and 99% confidence contour 
levels of the ISM density and SRCUT break frequency derived 
from the spectrum of region e. Lower panel: 68%, 90%, and 99% 
confidence contour levels of the ISM density and the ISM ion- 
ization parameter derived from the spectrum of region e. The red 
curves are isochrones at 100 yr and 250 yr after the interaction 
with the shock front. 
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Mg, and Si abundances to the values derived in region e, where 
the contribution of thermal emission is maximum (in the other 
regions, the nonthermal continuum is higher). In regions / and 
g, we found that O and Ne are instead significantly different from 
those in region e. For example, in region g we get;^ 2 = 588.7 
(507 d. o. f.) by letting O and Ne free and x 2 = 636.7 (509 
d. o. f.) with O and Ne tied to the values of region e. According 
to the F-test, the probability that the improvement of the fit is 
not significant is < 2.5 x 10~ 9 . The best-fit values of O and 
Ne in regions e, /, g, h are shown in Table 1 . In region h, the 
abundances were tied to those of region g, because they were all 
consistent within 0.5 cr. We then verified that the ejecta temper- 
ature and ionization parameter do not vary significantly across 
the regions. Therefore, we fixed them to their average values 
(kT ejec ta = 0.397 keV and t = 5.29 x 10 10 s cm 4 , respectively). 

In all the regions, we detected the ISM component. The best- 
fit parameters for all eight spectral regions are given in Table Q] 



While we did not observe any clear trend in the T e azimuthal 
profile, we found that the postshock density significantly in- 
creases in regions where particle acceleration is more efficient 
(Table Q] and Fig. 0J. Conside ring that the presho ck density is 
expected to be fairly uniform (Dubne r~ st al.ll2002h . the density 
profile in Fig.|4]shows that the shock compression ratio is higher 
near the nonthermal limbs. This result looks consistent with the 
predictions of the shock modification theory. 

It has been shown that where particle acceleration is effi- 
cient, models that do not account self-consistently for the ac- 
celeration process (as VPSHOCK) can provide biased results 
dPatnaude et al.l l2010h . We did not analyze spectra extracted 
from the nonthermal limbs, and we do not expect that parti- 
cle acceleration is at its maximum efficiency in regions a - h. 
However, some effects may be present in the regions closer to 
the nonthermal limbs, wher e "moderate" accelerat ion is at work 
(see Sect. @}. Nevertheless. iPatnaude et all ( 120101) showed that, 
in conditions of moderate particle acceleration, models that do 
not account for particle acceleration provide inaccurate results 
below 1 keV. They also showed that the main problems are vis- 
ible with high statistics (much higher than that of our observa- 
tions) at energies corresponding to the O lines, which are not 
well modelled (this issue may somehow affect the estimation of 
the ionization parameter). We can reasonably conclude that our 
estimates of the density of the ISM component, whose contribu- 
tion is higher at energies > 2 keV (see Fig. |2|, are reliable and 
that the density enhancement toward the nonthermal limbs is ro- 
bust. A more detailed comparison between the available spectral 
models and models including the effects of the acceleration pro- 
cess is beyond the scope of this work, but is required to address 
this issue. 

We point out that the values of the ISM density (derived by 
the ISM emission measure as explained in Sect. [2) do not change 
significantly by relaxing the assumption of pressure equilibrium 
between ISM and ejecta. Indeed, the increase in the ISM density 
toward the nonthermal limbs shown in Fig.|4]is also recovered by 
assuming a more realistic scenario, where the ejecta pressure is 
lower than the ISM pressure. Detailed 3 -D hydrodynamic mod- 
eling of SN 1006 dOrlando et alJl2012l) shows that the average 
ratio of the ejecta to the ISM pressure is P e j/PisM = Prei ~ 0.3. 
T herefore, as a further c heck, we derived the ISM densities as 
in lBocchino et alJ (1 19991) by assuming P re i = 0.3 and obtained 
the results shown in the upper panel of Fig. [5] The lower panel 
of the same figure shows that, even in the unrealistic scenario of 
P re \ = 2 (i. e. the ejecta pressure is twice the ISM pressure), the 
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Figure 4. Azimuthal profile of the postshock ISM density (see 
regions in Fig. 1). Error bars are shown at the 68% (crosses) and 
90% (shaded areas) confidence levels. The blue curves show the 
profiles predicted bv lVolk et ail (l2003h with SB/B=1 (dot) and 
0.8 (dashed), and derived by the radio flux (solid curve, see text) 
by assuming the relationship between t he injection ef f icienc y 
and shock compression ratio adopted by Ferra nd et al.l {2010). 
The green line indicates the expected trend in the case of no 
shock modification. 



ISM density near the nonthermal limbs is significantly higher 
than in the center of the thermal limb. 

We found that the filling factor of the ISM, f ISM , in all the 
spectral regions is higher than 95%. Hence, the volume fraction 
occupied by the ejecta is relatively small. This is in qualitative 
agreement with detailed magneto-hydrodynamic simulations of 
SN 1006 that show t he presence of narr ow fingers of ejecta close 
to the shock front (Orlan do et al.ll2012l) . 

Table 1 shows that the ejecta density is somehow correlated 
with the ISM density. This may suggest that effects of shock 
modification are present also at the reverse shock. However, 
while our estimates of the ISM density are robust, the ejecta den- 
sities may be less reliable, since they strongly depend on our es- 
timates of the filling factors. The small relative variations of /ism 
in regions a- h (/ism ~ 0.95 - 0.98), in fact, correspond to large 
relative variations of the ejecta filling factor (//« ,« ~ 0.02-0.05) 
that strongly affect the estimates of the ejecta density. For exam- 
ple, by assuming that the filling factor of the ejecta is the same 
in all the regions and is equal to its average < f e j ecta >= 0.031 
(and P re i =1), the ejecta density in region a and e are n e j ~ 0.52 
cirT 3 and n e j ~ 0.60 crrT 3 , respectively (to be compared with the 
values reported in Table [TJ. On the other hand, we verified that 
the estimates of the ISM density are not affected by this issue. 

Table 1 also shows the values of the cut-off frequency Vbreak, 
which increase s toward the nont hermal limbs, as expected. 
With respect to iMiceli et al.l (120091) (where larger regions were 
adopted), our values are slightly smaller. This effect is due to the 
prese nce of the ISM component (not included in IMiceli et al.l 
2009). Because of its extreme underionization (t < 7 x 10 8 
s cm -3 ), the ISM component is almost featureless, and in the 
archive observations where the statistics were not very high, it 
was possible to model its contribution by increasing the cut-off 
frequency of the SRCUT component. 

The values of the ejecta ionization parameter in all regions 
is compatible with ~ 5 x 10 10 s cm 4 . These values are a fac- 



5 



M. Miceli et al.: XMM-Newton evidence of shocked ISM in SN 1006: indications of hadronic acceleration 



Table 1. Best-fit parameters 



Parameter 

nn-i -J% 


Region a 

i nrv+0.4 


Region b 

i ni+U.OS 


Region c 
i 01 +0.0!) 


Region d 

1 tO+U.UK 


Region e 


Region / 


Region g 


Region h 



£r ISM (keV) 

T ISM (10 8 scm- 3 ) 

n ejeaa (10" 'cm -3 ) 

O 

Ne 

Mg 

Si 

v bmk (10 15 Hz) 
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All errors at 90% confidence level. kT e j eaa and T e s eaa are frozen to 0.397 keV and 5.29 X 10 s cm , respectively. 
The densities are estimated from the emission measure by calculating the volume of the emitting plasma and by 
assuming pressure equilibrium between ejecta and ISM. 
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Figure 5. Upper panel: ISM densities derived in regions a - h of 
Fig. 1 by assuming that ratio of the ejecta to the ISM pressure is 
P re i = 0.3. Errors at 90% confidence level. Lower panel: Same 
as upper panel for P re i = 2. 



tor ~ 2-3 higher than those reported in Mic eli etaD (|2009). 
We verified that these discrepancies are associated with the up- 
graded version of the VPSHOCK model (based on APED) that 
we are adopting here. Since the ejecta densities are ~ 0.6 cm -3 , 
we derive that the plasma was shocked t s i, oc k ~ 2500 yr ago, i. e., 
more than the age of SN 1006. This puzzling result relies on the 



assumption that the density remains constant over time. Lower 
values of t s h oc k can be obtained by considering that the ejecta 
density decreases during the evolution of the remnant (because 
of the expansion). A lower value of t s hock can a l so be obtained by 
assuming that the ejecta component consists of a "pure ejecta" 
plasma, where the number of free electrons is much higher than 
the number of positive ions (because heavy elements are much 
more abundant than hydrogen). We verified that by increasing 
the chemical abundances by a factor of 1000 in the ejecta com- 
ponent and by reducing the emission measure n e nnV by a similar 
factor (with respect to the values that can be derived from Table 
1) we still obtain a good fit in region e. This indicates that a pure 
ejecta scenario is plausible. In this case, the ejecta emission mea- 
sure (~ 1000n e n//V) would be almost the same as that reported 
in Table 1 , but the electron density would be higher by a factor 
Vz (where Z is the average number of free electrons per ion) 
and the time elapsed after the shock impact would be reduced 
by the same factor. We will investigate in detail this scenario in 
a forthcoming paper (Miceli et al., in preparation); we here only 
note that, even within this pure ejecta model, the normalization 
of the ISM component does not vary significantly. 



4. Discussion and conclusions 

We have analyzed a set of deep XMM-Newton observations of 
SN 1006 that have allowed us to address important issues. 

We have detected X-ray emission from the shocked ISM in 
the southeastern rim of SN 1006. It was not possible to explore 
the conditions of the shocked ISM in the nonthermal limbs be- 
cause the higher contribution of the synchrotron radiation did not 
allow us to obtain tight constraints. Nevertheless, in regions a—h, 
the high statistics of the new observations showed that the addi- 
tion of a thermal, underionized component improves the quality 
of our spectral modeling. A further indication that the associa- 
tion of this additional component with the ISM is sound is that its 
contribution to the nitrogen line emission is significant. Nitrogen 
line emission cannot be associated with ejecta in remnants of 
Type la supernovae and should originate in the shocked ISM. As 
shown in Fig. [2] the contribution of the ejecta component dom- 
inates below 0.5 keV (i. e., the band of the N lines). However, 
we found that the fit of our model to the spectra does not change 
significantly if we set the nitrogen abundance in the ejecta com- 
ponent to zero, while by putting N=0 in the ISM component, we 
obtained significant residuals around the nitrogen line. This in- 
dicates that in this band the contribution of the ISM component 
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to the nitrogen emission lines is, as expected, stronger then the 
contribution of the ejecta component to the same line. 

Figures [4] and [5] clearly indicate that, under the assumption 
of constant preshock density, the shock compression ratio in- 
creases (up to ~ 5.5-6) in regions where nonthermal emission 
is stronger. This is qualitatively consistent with the expectations 
from modified shock models. In particular, such a high shock 
compression ratio indicates that the fr actional cosmic- ray pres- 
sure in the postshock region is ~ 30% (IVinketalJl2010h near the 
nonthermal limbs. 

Though a detailed theoretical modeling of the expected az- 
imuthal profile of the shock compression ratio in SN 1006 is not 
available in the literature, we compared the observed ISM den- 
sity profile wit h that pred i cted b y the shock modification model 
provided by IVolk et all d2003l) and shown as dashed/dotted 
blue curves in Fig. [4] The solid blue curve was instead de- 
rived by assuming that the injection efficiencjQ, rj, is propor- 
tional to the radio flux divided by B 3/2 (and has its minimum 
at rj = 5 x 10 s in region e ) , and adopting the magnetic field 
model MF2 of iPetruk et all (2009). The green line shows the 
constant compression-ratio scenario, which corresponds to no 
shock modification. For all the curves, the relationship between 
the injection efficiency and the shock compressio n factor has 
been obtained by following iFerrand et alj d2010h . The shock 
modification models predict an azimuthal trend similar to the 
measured one, though they cannot fit all the points. In particular, 
the models shown in Fig. [4] slightly underpredict the observed 
density variations. However, we point out that a complete, self- 
consistent model of hadron acceleration at oblique shocks is not 
available yet and that current models do not provide tight theo- 
retical constraints, though they allow us to get qualitative indi- 
cations. Moreover, the models available do not include a depen- 
dence of the compression ratio on the maximum energy of the 
accelerated hadrons which may be higher near the nonthermal 
limbs. The cut-off frequency of the synchrotron emission from 
ultrarelativistic electrons is, in fact, much hig her in the limbs 
than i n the center and elsewhere in the rim dRothenflu g et alj 
2004: lMicelietalJl2 009). This indicates that electrons are accel- 
erated to higher energies in the limbs and suggests that protons 
may also be accelerated to higher energies. 

In conclusion, we have detected the contribution of the 
shocked ISM to the X-ray emission in the southestern rim of 
SN 1006. We also found an azimuthal trend of the ISM post- 
shock density that suggests that cosmic ray acceleration is at 
work at the shock front of SN 1006 and modifies the shock com- 
pression ratio. The results we have obtained in the southeastern 
rim of SN 1006 therefore indicate the presence of shock modi- 
fication by particle acceleration. A definitive proof will require 
more data at higher spatial resolution. These data are necessary 
to extract spectra from narrow regions between the shock front 
and the contact discontinuity in order to better isolate the contri- 
bution of the shocked ISM and to confirm the effect that we have 
detected in this work. 
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Appendix A: Volume calculation 

We here describe the procedure we developed to calculate the 
volume of the emitting plasma in the spectral regions shown in 
Fig. 1. 




Figure A.l. Schematic representation of SN 1006 (sphere with 
radius R s heii) and of the spectral regions shown in Fig. 1 (yellow 
area) 



Appendix B: Analysis of the pn spectra 

To further test the robustness of the density profile shown in Fig. 
|U we also fitted the pn spectra only in the 0.6 - 7 keV band. 
Alhough in this case the best-fit parameters are more sensitive to 
the background subtraction procedure (specially the ISM tem- 
perature) and error bars are larger, we still obtained an ISM den- 
sity profile consistent with that of Fig. [4] Figure IBTI shows the 
density profile obtained by fitting the pn spectra only and by 
adopting method 1 (see Sect. |2| for the background subtraction. 
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Let us assume that the SN 1006 shell is a sphere with ra- 
dius R s heii and let us consider a cylinder with radius Rj nt an d 
height 2R s heu- In agreement with Cassam- Chenai' et"aT] ([2008), 
we consider as the center of the shell the point with coordi- 
nates a jmo = 15 /, 02'"56.8 s ,<Jy20oo) = -41°56'56.6") (see Fig. 
1). We also assume a distance of 2.2 kpc. Let R s heiu Rim, and 9 
be the outer radius, inner radius, and opening angle of our spec- 
tral regions, respectively. Figure lATI shows the projection of the 
sphere and the cylinder in the plane of the sky (left panel) and 
in the plane identified by the line of sight and the east direction 
(right panel). Since the volume of the spherical cap with height 
h = RsheU-(R 2 shd -R 2 int ) XI2 is V cap = 7i h / '6(3 R 2 int +h 2 ), the volume 
of the intersection between the cylinder and the sphere is 

Vcs = 7iR 2 nt x 2R sheU - 2(nR 2 nt h - V cap ) (A. 1) 

The volume of the portion of the sphere outside of the cylin- 
der is Vo = 4/37rR^ ;ifj// - V c - S and the volume of our spectral 
region is V reg = 8/2jtVo which depends only on the outer radius, 
inner radius, and opening angle of the region. 

For each region, we then measured R s heii (SN 1006 is not 
a perfect circle in the sky and small azimuthal variations are 
present), /?,„,, and 8 for each region, thus obtaining V reg - a 
23 x 10 56 cm 3 , V reg - b = 2.7 x 10 56 cm 3 , V reg - C = 3.4 x 10 56 
cm 3 , V reg - d = 3.0 x 10 56 cm 3 , and V reg - e = 3.2 x 10 56 cm 3 . 



Figure B.l. ISM densities derived in regions a - h of Fig. 1 by 
fitting the pn spectra only in the 0.6 - 7 keV band. Errors at 90% 
confidence level. 
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